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Steady thermo-solutocapillary convection in a rectangular cavity with deformable free surface under microgravity 


condition is numerically studied, where level set method is employed to capture the free surface deformation. 


Both the temperature and solute concentration gradients are applied horizontally. The computational results show 


that, as the thermal to solutal Marangoni number ratio varies between -10 and -1 (namely, -10 < Rg < -1), the flow 


field exists one anti-clockwise rotating convective cell driven by thermocapillary convection, and the free surface 


bulges out near the left end wall and bulges in near the right end wall. As -1 < Rg < 0, the flow field exists one 


clockwise rotating convective cell driven by solutocapillary convection, and the free surface bulges out near the 


right end wall and bulges in near the left end wall. As Rg = 


-1, the flow field consists of one clockwise and one 


anti-clockwise rotating convective cells, and the free surface bulges in at the central point and bulges out near the 


left and right end walls. 
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Introduction 


Thermo-solutocapillary convection is driven by sur- 
face tension gradient due to the variation of temperature 
and solute concentration at free surface. This phenome- 
non widely exists in crystal growth process, welting pro- 
cess, and many other industrial processes. The most 
widely studied model associated with the thermo-soluto- 
capillary convection is a rectangular cavity with temper- 
ature and solute concentration differences between two 
end walls. At 1986, Bergman [1] numerically simulated 
the double-diffusive Marangoni convection in a rectan- 
gular cavity with horizontal temperature and concentra- 
tion gradients, and he found that even though the oppos- 
ing thermal and solutal capillary forces are of equal mag- 
nitude, convection occurs when thermal Marangoni nu- 
mber exceeds a critical number. Chen et al [2] studied the 
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onset of steady and oscillatory flow of double-diffusive 
Marangoni convection when the resultant thermal and 
solutal Marangoni effects are equal and opposing, and the 
effects of the Lewis number, the Prandtl number and the 
aspect ratio on stability were analyzed. Similarly, Li et al. 
[3] further investigated the transitions to chaos in double- 
diffusive Marangoni convection in a rectangular cavity 
with horizontal temperature and concentration gradients. 
Zhan et al. [4] numerically investigated three-dimen- 
sional double-diffusive Marangoni convection in a cubic 
cavity with both the temperature and solute concentration 
gradients were applied horizontally, and the effect of 
surface-tension Reynolds number, surface tension ratio 
and Lewis number on the flow regime was discussed. 
Chen and Chan [5] studied the stability of double-diffu- 
sive Marangoni convection in a horizontal fluid layer 
driven by the combined effects of buoyancy and surface 
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Nomenclature 
C solute concentration, wt% 
C 
a thermal Capillary number, Ca; = Fral 
O0 
Ca . A 
c solute Capillary number, Cac = ta 
O0 
solute diffusivity (m*-s”) 
k curvature 
L characteristic length (m) 
Le 
Lewis number, Le = = 
l 
Mac 
; ACL 
solutal Marangoni number, Maç = Ae 
MQ) 
Mar 
: ATL 
thermal Marangoni number, Map = Yr 
Ma] 
n unit vector normal to the interface 
Pressure (Pa) 
Pr Q) 
Prandtl number, Pr = — 
vı 
Re UL 
Reynolds number, Re = BES 
Hı 
R, ratio of thermal to solutal Maragoni number, 
_ Yr AT 
7  ¥¢+AC 
t time (s) 
T temperature (K) 
u velocity vector (m.s”) 


tension. Recently, Li et al. [6] investigate the characteris- 
tics of the coupled thermal and solutal capillary convec- 
tion with the radial temperature and solute concentration 
gradients in a shallow annular pool by asymptotical 
analysis. Furthermore, Li et al. [7] simulated the transi- 
tion characteristics of the thermal-solutal capillary con- 
vection numerically in an annular pool. However, all 
these studies were dedicated to the situation of infinite 
small capillary number, and the deformation of free sur- 
face or interface was neglected in order to avoid the si- 
mulation complexity associated with the unknown free 
surface or interface shapes. 

It is well known that the deformation of free surface or 
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U characteristic velocity, 
U = Max(y,AT}|ycAC|)/ s (m's") 

x,y coordinate system (m) 

Greek symbols 

a thermal diffusivity(m7.s”) 

ġ level set function 

Ve surface tension temperature coefficient 
(N-m'T”) 

ie surface tension solute concentration coeffi- 
cient (N-m'-wt%"') 

A thermal conductivity (W-m'': K’) 

u dynamic viscosity (kg-s'-m'') 

v kinematic viscosity (m*-s"') 

p density (kg'm”°) 

o surface tension, 
o=09 + q(T -T)+rc(C- Co) Nm’) 

Oo initial surface tension at T = Tọ and C = Cy 

Subscripts 

g gas 

l liquid 


interface has an important influence on the pure thermo- 
capillary convection, in particular for the onset of the 
oscillatory flow [8]. Microgravity experiment by Koster 
[9] indicated that, a detailed investigation of thermoca- 
pillary convection in multilayered fluid systems needed 
to take account for finite interface deformations in con- 
junction with a suitable dynamic contact line condition. 
Moreover, the deformation effects of free surface or in- 
terface in the numerical simulation of pure thermocapil- 
lary convection was considered by many investigators 
[10-14]. 

The authors first considered the dynamic deformation 
of free surface in the numerical simulation of thermo- 
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solutocapillary convection in our previous paper [15], 
and in the present work we further investigate the influ- 
ence of thermal and solutal Marangoni effects on free 
surface deformation in an open rectangular cavity. 


Physical and Mathematical Models 


The two-dimensional model is an opened liquid-filled 
rectangular cavity, as shown in Figure 1, and the length 
and height are L and H, respectively. The cavity is filled 
with a binary fluid, and the gas-liquid interface on the top 
is a deformable free surface. Different temperatures and 
concentrations are applied at the left (Tı, Cı ) and right 
(T2, C2) end walls, where Tı < T2 and C; < C2. The non- 
slip boundary condition is adopted for the end and bot- 
tom walls, and the surface tension force is applied on the 
free surface. All the thermal properties are assumed to be 
constant except for the surface tension, which is allowed 
to vary linearly with the liquid temperature and solute 
concentration. Thus the surface tension can be formu- 
lated as, 


O=0o4 yr(T To) H Yc(C Co) 
oO (To, Co), Y= —(00/6T).. à 


Thermophysical properties of the 


Where op = 


yc = —(00/0C),, . 
fluid are estimated at the reference temperature Ty and 
concentration Co, which are set to be equal to T, and Cj, 
respectively. Moreover, only the liquid surface tension 
increases with concentration and decreases with temper- 
ature are considered. The fluid flow is assumed to be 
laminar, and buoyancy effects are neglected in the 
present study. 


gas 
Ti, ol To, 
0 


Fig. 1 Physical model 


The non-dimensional variables of velocity, tempera- 
ture, concentration, coordinates, time and pressure are 


defined as V a os = C= a , (X,Y) 
U T, -T, C-C, 
a, T -2 P = pl(pU"), respectively. Here, U = 


Max(|yrAT 


YCAC I) / 44 is characteristic velocity. Then, 


5) 


the non-dimensional Navier-Stokes equations governing 
the thermo-solutocapillary convection are expressed as 
follows: 

V-V=0 (1) 
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OV 1 1 2 
—+V-(VWV)=-—VP+——V-(LVV 
ape) BO eRe (avv) 
1 r 1l VC 
V-(aVV ——| VO 5 ô 2 
TV +h Epo @ 


PRe\ Ca Ro a) 
® .v.-ve-—1—Vv (ive) (3) 
OT pCMa 
A sv.v = _v.(Bvc") (4) 
OT PrLeRe 


where f= u/1,,P=p/p,, &4=4/4, D=D/D, 
and C p =C,/C,) are the dimensionless viscosity, den- 
sity, thermal conductivity, solute diffusivity, and specific 
heat, respectively. 5 is the smeared-out Dirac delta func- 
tion, and V,=(I-nn)-V is free surface gradient oper- 


ator. The dimensionless parameters are defined as fol- 
_ YrATL 


lows: Thermal Marangoni number Ma, = , So- 
MQ) 
; ACL : 
lutal Marangoni number Maç sAn , Marangoni 
Hy) 


number Ma=Max(| Mar, |, | Mac |), Reynolds number 


kee Capillary number Ca = Max(|Ca;|, | Cac |), 
ri 

the ratio of thermal Marangoni number to solutal Ma- 

ma , Prandtl number Pr = 2 and 

Yc: AC Vy 


rangoni R, = 


Lewis number Le =£. The subscripts of g and 1 note 


i 
the gas and liquid, respectively. The continuum surface 
force (CSF) model [16] is used to reformulate the surface 
tension as a volume force, and the fourth and fifth terms 
at the right hand of the momentum equation are the tan- 
gential and normal surface tension, respectively. 

The level set method is employed to capture the free 
surface implicitly by introducing a smooth level set func- 
tion, with the zero level set as the free surface, positive 
value underside the free surface, and negative value up- 
side the free surface. Consider the following free surface 
evolution equation: 

Og 
—+V-Vd=0 (5) 
OT 

which will evolve the zero level of ø= 0 exactly as the 
actual free surface moves. The corresponding physical 
variants can be expressed as 


b=n,+(1-n,)H(9) 


(6) 
fi=n, +(1-m,)H(9) 
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Where 77, =P,/Pp Nu = Heki, H() is the smeared-out 
Heaviside function defined by 


0 
Jlj¢ ta (@) #8 
a= ares daag 
1 o>é 


where ¢é is a tunable parameter that determines the size 
of the bandwidth of numerical smearing. A typical good 
value is é= 1.5AY. 

To keep the level set function as a distance function 
from the front, an approach based on solving the hyper- 
bolic partial differential equation has been presented in 
reference (Sussman et al. [17]). The reinitialization equa- 
tion 1S 


dy =sign(%)(1-|V9]). 9.0) =A) ® 


Where \Vdl= \ de + o , and the sign function 
g 


Sela T 


Moreover, in order to make the total mass completely 
satisfy the mass conservation in time, the mass conserv- 
ing procedure proposed by Liang et al. [18] is used. 

Boundary and Initial Conditions 


e=-0,c-0,v=-0, E atx-0 (9-a) 
ax 
og 

@=1, C=1, V=0, 0 at X=1 (9-b) 
ax 

a0 3h, 2 sy: ©! =p (9-c) 

aY aY Y 

@=0, C=0, V=0, ø(X,A)=0,P=0att=0 (9-d) 


In this paper, we consider the contact conditions of 
interface with the end walls is contact points fixed, and 
the gas-liquid interface is a thermally insulated interface. 


Computational method 


In this article we use the RKCN projection method 
presented by Ni et al. [19] to solve the control equations. 
The Crank-Nicholson implicit technique is employed to 
update the diffusion term, and the low storage three-stage 
Runge-Kutta technique is employed to update the con- 
vective term. This method has second-order temporal 
accuracy for variable-density unsteady incompressible 
flows. The diffusion term in momentum equation is dis- 
cretized using standard central difference schemes, while 
the convective term is discretized by high-order compact 
schemes. The convective term in the level set equation 
using the third-order ENO scheme. The time step for all 
the simulations is 1 x 10°. This detailed description of 
the computational method can refer to our previous paper 
[20]. The numerical code is validated by comparing 
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computed surface deformation of pure thermocapillary 
convection with those from Sim [8] in Figure 2. The pa- 
rameters for the simulation are Pr = 1, Ar = 1, Ca = 0.05 
and Re = 1. The numerical results at Bi = 0 are in good 
agreement with Sim’s results. 

The two-dimensional rectangular cavity is discretized 
with a uniform mesh, and this mesh is determined by 
refining the mesh size until convergence and constancy 
of flow velocity, temperature, concentration and free sur- 
face deformation is achieved. In this paper the grid 400 x 
x 200 y is adopted, which meets the demands of conver- 
gence. 


Sim's results 
Present results 


0 
-0.0005 0 0.0005 


Fig. 2 Comparison of free surface deformation with Sim’s 
results 


Results and Discussions 


In the open double-diffusive fluid cavity, both the 
temperature and solute concentration differences between 
the left and right end walls produce a surface tension 
gradient at free surface, so the thermocapillary and solu- 
tocapillary convection in the liquid layer can be induced. 
The relative importance of the thermocapillary and solu- 
tocapillary convection for the flow depends on the ther- 
mal to solutal Marangoni number ratio R,, therefore, 
different convective cell structures and free surface de- 
formation modes under different thermal to solutal Ma- 
rangoni number R, can occur. The numerical simulations 
are carried out for Prandtl number Pr = 1, Capillary 
number Ca = 0.1, Marangoni number Ma = 10, Lewis 
number Le = 10, and thermal to solutal Marangoni num- 
ber ratio -0.1<R,<-10. 

At first, the case of thermal to solutal Marangoni 
number ratio R, = -10 is considered, which means the 
thermocapillary convection dominates the fluid flow. 
Figure 3 shows the streamline, temperature and solute 
concentration field distribution. It can be seen that the 
flow structure has an anti-clockwise rotating convective 
cell and occupies the whole cavity. Under the action of 
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dominant thermocapillary forces at free surface, the free 
surface bulges out near the left wall and bulges in near 
the right wall. At the same time, the temperature and so- 
lute concentration contours manifest as gradient distribu- 
tion along the horizontal direction, this is due to the con- 
vective intensity is weaker at Ma = 10 and the flow is a 
steady flow when the thermo-solutocapillary convection 
fully developed. Moreover, the distortion intensity of 
solute concentration contours is larger than that of the 
temperature contours, that is because the solute diffusiv- 
ity is smaller than thermal diffusivity as Le = 10. 

Due to the thermo-solutocapillary convection is driven 
by surface tension gradient, so the mechanism of flow 
structure and the free surface deformation mode can be 
interpreted by the surface tension distribution. Figure 4 
gives the variation of non-dimensional surface tension 
along horizontal direction at Ro = -10. In the figure, the 
non-dimensional surface tension is monotone decreasing 
with the increase of coordinate X. It is well known that 


(c) 


Fig. 3 Streamline, temperature, and solute concentration dis- 
tribution at R, = -10 ((a)streamlines, (b)temperature, 
and (c)solute concentration distribution) 
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0.96 
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Fig. 4 Variation of non-dimensional surface tension along 
horizontal direction at R,=-10 


the fluid flow driven by surface tension at free surface is 
from the smaller surface tension spot to the larger surface 
tension spot. Therefore, under the action of the monotone 
decreasing surface tension the flow at free surface is from 
the right end wall to the left wall, then leads to the fluid 
accumulates in the vicinity of the left end wall and the 
free surface bulges out. 

Figure 5 gives free surface deformation under differ- 
ent thermal to solutal Marangoni number ratios as -10 < 
R,< -1. These free surfaces bulge out near the left end 
wall, and bulge in near the right end wall. The free sur- 
faces is almost asymmetric about the point (X = 0.5, Y = 
0). With the increase of the thermal to solutal Marangoni 
number ratio the deformation intensity of free surface is 
decreased. The surface deformation is O(10°), and its 
maximum value is 5.4 x 10° with Ca = 0.1 Ma = 10 and 
Le = 10. Figure 6 shows horizontal velocity distribution 
at free surface under different thermal to solutal Maran- 
goni number ratios as -10 < Ro < -1. In the figure, the 
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Fig. 6 Horizontal velocity distribution at free surface as -10 < 
R,<l 
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horizontal velocity at the free surface is minus, and with 
the thermal to solutal Marangoni number ratio increasing 
the free surface velocity is decreased. That is because the 
fluid flow is driven by thermocapillary convection, and 
thermocapillary convection intensity is strengthened as 
R, decreases. 

Next, the case of thermal to solutal Marangoni number 
ratio R, = -0.1 is considered, which means the solutoca- 
pillary convection dominates the flow in the cavity. The 
streamlines, temperature and solute concentration con- 
tours distribution is shown in Figure 7. It can be seen that, 
the flow field has a clockwise rotating convective cell, 
this is reverse to that of the case R, = -10, this is due to 
the driving direction of solutocapillary forces opposite to 
that of thermocapillary forces. Moreover, the free surface 
bulges out near the right end wall and bulges in near the 
left wall. The temperature and solute concentration con- 
tours are gradient distribution along the horizontal direc- 
tion, and the distortion intensity of solute concentration 
contours is larger than that of the temperature contours. 

Figure 8 gives the variation of non-dimensional sur- 
face tension along horizontal direction as R,= -0.1. In 
the figure, the non-dimensional surface tension is mono- 
tone increasing with coordinate X. Therefore, under the 
action of the monotone increasing surface tension the 
flow at free surface is from the left end wall to the right 
end wall, so leads to the fluid accumulates near the right 
end wall and the free surface bulges out. 

Figure 9 gives free surface deformation under differ- 
ent thermal to solutal Marangoni number ratios as -1 < 
R,< -0.1. These free surfaces bulge out near the right end 
wall, and bulge in near the left wall. The free surfaces is 
almost asymmetric about the point (X = 0.5, Y = 0). The 
surface deformation is O(10°), and its maximum value is 
5.0 x 10° with Ca = 0.1 Ma = 10 and Le = 10. Moreover, 
with the increase of the thermal to solutal Marangoni 
number ratio the deformation intensity of free surface is 
increased. 


(c) 


Fig. 7 Streamline, temperature and solute concentration dis- 
tribution as R,= -0.1 ((a) streamlines, (b) temperature, 
and (c) solute concentration distribution) 
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Fig. 8 Variation of non-dimensional surface tension along 
horizontal direction as R, = -0.1 
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Fig.9 Free surface deformation as -1 < R,< -0.1 


Figure 10 shows horizontal velocity distribution at 
free surface under different thermal to solutal Marangoni 
number ratios as -1 < R,< -0.1. The horizontal velocity 
at the free surface is positive, and with the thermal to 
solutal Marangoni number ratio increasing the free sur- 
face velocity is increased. That is because the fluid flow 
is driven by solutocapillary convection, and its convec- 
tive intensity is strengthened as R, increases. 

Finally, the case of thermal to solutal Marangoni 
number ratio Rọ = -1 is considered, this means that the 
opposing thermocapillary and solutocapillary forces are 
of equal magnitude at free surface as the flow fully de- 
veloped, then the whole fluid system manifests as a bal- 
ance system. Figure 11 gives the flow, temperature, and 
solute concentration fields distributions as R, = -1. It can 
be seen that the flow field consists of two convective 
cells, where the left cell is anti-clockwise rotating and the 
right cell is clockwise rotating. That is because both 
convective cells are driven by thermocapillary forces and 
solutocapillary forces, respectively. Both the temperature 
and solute concentration contours are strictly gradient 
distribution along the horizontal direction due to the 
weaker convective intensity. 
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Fig. 11 Streamline, temperature and solute concentration dis- 
tribution at R, = -1 ((a)streamlines, (b)temperature, 
and (c)solute concentration distribution) 
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Fig. 12 Variation of surface tension along horizontal direction 
as R= -1 
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Fig. 13 Free surface deformation as R= -1 
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Fig.14 Horizontal velocity distribution at free surface as R, = -1 
Conclusions 


In this paper, thermo-solutocapillary convection in a 
rectangular cavity with dynamic deformable free surface 
under microgravity condition is numerically simulated, 
and the effect of thermal to solutal Marangoni number 
ratio on the flow and free surface deformation is dis- 
cussed. The following conclusions were drawn: 

1) As the thermal to solutal Marangoni number ratio 
varies between -10 and -1, namely -10 < R,< -1, the flow 
field exists one anti-clockwise rotating convective cell, 
which is driven by thermocapillary convection, and the 
free surface bulges out near the left end wall and bulges 
in near the right end wall. 

2) As the thermal to solutal Marangoni number ratio 
varies between -1 and 0, namely -1 < R, < 0, the flow 
field exists one clockwise rotating convective cell, which 
is driven by solutocapillary convection, and the free sur- 
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face bulges out near the right end wall and bulges in near 
the left end wall. 

3) As the thermal to solutal Marangoni number ratio 
equals to -1, namely R, = -1, the flow field consists of 
one clockwise and one anti-clockwise rotating convective 
cells, they are driven by solutocapillary and thermocapil- 
lary convection, respectively. The free surface bulges in 
at the central point and bulges out near the left and right 
end walls. 
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